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Hundredths of kilometers above us there is a stormy 
World sea: the ionosphere, through which our EME signals 
have to go trough twice. The interactions are of various 
Co nfe re nce types: there are significant effects on signal amplitude 
(QSB), and wave polarization rotations (Faraday effect). 
Starting from data obtained from MAP65 decodes, we 


have made an analysis of QSB showing its dependence 
IONOSPHERIC 


not on attenuation but on focusing or defocusing 


INTERACTIONS effects due to ionospheric waves. Then we have 
WITH EM F SIGNALS defined the algorithms for calculating Faraday rotation 


over the total moon pass, and compared them to actual 
IK1 UWL-IK3XTV MAP 65 decodes. We have found also how Faraday 
typically behaves as a function of correspondent station 
orientation. All this analysis is focused mainly on 144 
MHz EME (our experience), but a good part is applicable 
to all bands. 
Note: this document contains the complete 
presentation, plus some additional detail information. 


IONOSPHERIC INTERACTIONS WITH EME SIGNALS 


Giorgio Marchi, IK1UWL — Flavio Egano, IK3XTV - EME 2014 Conference Date: August 2014 


Study start 
Ham EME communication through the ionosphere has evolved over the years, from CW to the first digital 


mode JT44, then the ubiquitous JT65, and today MAP65 (thanks Joe Taylor). 

Latest software arrival for EME communication, for a station equipped with cross yagis and suitable 
hardware, MAP65 gives JT65 type message decodes over an 80 kHz wide band, adding info on level and 
polarization. But MAP65 is not only a communications method, it is also a scientific instrument which, 
besides message decodes, gives us two important data: signal level and polarization angle. 

An EME pile-up 

Our interest for these two additional parameters arose from the beginning. Giorgio IK1 UWL had just installed 
MAP65 and used it to monitor OX3LX (and qso him). He saved 46’ of data 

UTC Date: 2012-ago-03 00:00 


QRG DF DT Pol dB UTC 

144.143 -129 0-2 0 1.7 63 4-22 0008 CQ OX3LX HP15 110 8 
144.143 -138 3-1 0 1.7 64 3 -200010 CQ OX3LX HP15 11015 
144.143 -144 0-2-1 1.7 48 3-22 0012 CQ OX3LX HP15 1105 
144.143 -153 0-2-1 1.9 43 5-22 0014 CQ OX3LX HP15 110 3 
144.143 -161 1-1 0 1.5 31 1 -210016 CQ OX8LX HP15 110 4 
144.143 -170 0-1 0 1.7 15 2-210018 F6HVK OX3LX HP15000 10 5 
144.143 -176 0 00 3.6 0 4-180020 RRR 000 

144.143 -185 10 0 1.7 18 5-200022 RK3FG OX8LX HP15000 1 017 
144.143 -191 0 00 1.0 0 4-190026 RRR 000 

144.143 -199 1-1 0 1.7 10 4-19 0028 CQ OX3LX HP15 110 3 
144.143 -205 0-2 0 1.5 4 4-180030 ISMEK OX3LX HP15 000 1 016 
144.143 -214 0 0 0 3.6 0 3-18 0032 RRR 000 

144.143 -217-1 0-1 2.1 0 4-180034 IZ3KGJ OX3LX HP15 000 1 018 
144.143 -226 0 00 1.0 0 4-23 0038 RRR 000 

144.143 -229 -1-1-1 1.6 -4 4-190040 CQ OX3LX HP15 11015 
144.143 -232 -2-1-1 1.8 -7 5-20 0042 CQ OX3LX HP15 11012 
144.143 -238 0-1-1 1.8 -9 4-180044 CQ OX3LX HP15 11010 
144.143 -243 3-1 1 1.8-13 3-200046 IKIUWL OX8LX HP15 O00 10 7 
144.143 -246 000 1.0 0 4-180048 RRR 000 

144.143 -252 1-1 1 1.8 23 6 -220050 CQ OX3LX HP15 1105 
144.143 -255 -1 -1-1 1.8 31 5 -22 0052 CQ OX3LX HP15 110 4 
144.143 -261 2 0-1 1.8 26 3-22 0054 CQ OX3LX HP15 1105 


and put them in diagram 
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Besides being glad for the rare dx, the amount of variation of signal level and polarization angle in so short a 


time was surprising. Giorgio IK1UWL sent this data to Flavio IK3XTV, and we both decided to investigate on 


the causes of this. 


Height Gan) 


The ionosphere 
- Partially ionized gas layer between ~50 and 


~1000 km height and permeated by Earth’s 
magnetic field is a turbulent ocean, roughened 
by high speed winds. 

- Free electron density is variable in space and 
time. 

- Their density (number per m3) determines various 
effects: 

Slowdown, Attenuation, Deviation, Rotation 


The ionosphere, space weather 
Weather near the surface involves winds that propagate like waves across the planet. 


Weather in the upper atmosphere and in space is characterised by very strong winds, so has wavelike 
characteristics as well. 

Space weather has two drivers: one is solar wind and radiation, the other is the dynamic atmosphere of 
Earth. 

In the atmosphere travel waves very low in frequency (wavelength in the range of hundreds or thousands of 
kilometres) that carry energy and momentum also upwards (buoyancy waves). 

Thermal effects cause tides with vertical movements and transfer of energy through shear winds. At the 
altitude at which sun radiation causes ionization there are winds with speed ranging in the hundreds of 


meters per second Being shear winds the ionosphere is roughened by shorter waves and vortices. 


= Jet stream in the troposphere, always from west 
= Jet stream in the mesosphere (reverses seasonally) 
= zonal currents in the stratosphere with seasonal inversion 


Eastward Winds — wetrom west 


Westward Winds _ eE=trom east 


0° 
‘ Equator 
Summer hemisphere Winter hemisphere 


Better position for Es 


Courtesy: Istituto di Scienze dell'Atmosfera e del Clima — CNR, Italia. 


These winds, which are active most days for part of the day, cause Travelling lonospheric Disturbancies 


(TIDs) which result in fluctuations of electron density and formation of amasses. 


lonospheric waves 
The ionospheric winds cause continuous undulations and waves (Traveling lonospheric Disturbances, TIDs) 


resulting in fluctuations of electron density. 
Horizontal Period Horizontal 
Class wavelenght — phase velocities 


LSTIDs | 51000 Km 0.5..3h | 300.1000 mis 
Large scale 


MSTIDs | 490.1000 Km |12 min...th | 100..300 m/s 
Medium scale 


SSTIDs <100 Km A few 
Small Scale minutes 


<200 m/s 


Source: INGV Istituto Nazionale di Geofisica e Vulcanologia - Italy 


This turbulence and waves in the ionosphere have a lens effect on our signals, 
focusing and defocusing them. 
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(Image source:Research and Technology Organisation. North Atlantic Treaty Organisation. Characterising 
the lonosphere. Author: G. Wyman (January 2009)) 


Rapid scintillations caused by ionospheric irregularities (ssTIDs, periods of some minutes) 


defocusing focusing 
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Slower scintillations due to msTIDs 
(observed at mid latitudes every day, wavelength 300 km, wind 100 m/s=360 km/h, period 50 minutes) 


QSB 
QSB band dependence (ionospheric refraction is proportional to inverse square of frequency) 
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Fig. 8. Detrended VHF, UHF, and L-band strip charts from a pass recorded at Poker Flat on 29 May 1976. 


Regions dominating effect 
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Signal level fluctuations 


A -— Static ionosphere 


Among the first notions of ionospheric behaviour we got, was that absorption of signal energy was in the 


range of 0,1 to 0,5 dB at 144 MHz. 


2 
i 
= 
é 
: 


lonospheric Absorption vs Frequency 


Typical daytime conditions 
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Image source: Darrel Emerson National Radio Astronomy Observatory - Tucson, AZ - USA 


B — Dynamic ionosphere 
Whilst the observed medium period fluctuations of signal level are in the order of 3-6 dB in day conditions 
(attributable to ssTIDs) 
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And of the order of 2 dB in night conditions 
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Additionally we noted long term fluctuations (attributable to msTIDs): 
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These cannot be attributed to absorption because a thirty-fold density change would be needed. 


So we realized that turbulence and waves in the ionosphere have a 
lens effect on our signals, focusing and defocusing them. 


Faraday effect, 
Faraday rotation is what happens to the polarization plane of a wave going through a ionized medium 


intersected by a magnetic field. The formula is very simple: 

@=k*B* TEC/f? (rad), 

with 

B = Geomagnetic field component in wave’s direction of travel (i.e. towards or from the Moon) 
TEC = Total Electron Content of the path 


f = wave frequency 


Effects tied to frequency 
Angle of rotation is inversely proportional to the square of frequency. 


Since our bands have frequencies increasing by a factor of three, the rotation becomes one ninth from a 


band to the next. 


With the same conditions for B and TEC: 


Band Rotation (9 

50 MHz 90 360 2.25 turns 
144 MHz 10 40 90 
432 MHz 1,1 4,5 10 
1296 MHz 0,1 0,5 1,1 


Evidently, Faraday is a concern mainly in VHF 

Microwavers are concerned only by Spatial Offset. 

Clarification on spatial offset: 

Polar polarization is the angle between an antenna and earth’s polar axis. 

Spatial offset between two stations is simply the difference between the polar polarizations of the two 
stations. This offset adds algebraically to Faraday rotation. 


Geomagnetic field 
We found a good source for F (total field) in the web site of the British Geological Survey. 


Introducing Lat, Long, Median Height of the ionosphere, and Date, one obtains: 


Total field F (nTesla) 
Inclination! (9 


Declination D (9 
Magnetic latitude 


Down 


> Geodetic Coordinates 


Latitude: [43.82 | degrees |O minutes fo | seconds 


Longitude: |7 83 degrees |O minutes |O | seconds 


Altitude: km above MSL. 


Date 
Date: (2013-01-20 | (Y¥YYY.YY OR YYY¥-MM-DD) 


Show result on map 


[mappa _]_ Satellite 
ae 


Latitude | Longitude | Altitude Date 
43.820 7.830 270.00 [2013.05 


Comp D 1 x | ¥ 
MF | 0.973 | 59.408 | 20904 | 355 | 20007 | s5s63 | 41084 
sy 65 0.3 137 [sos | 144 | 176 [ 225 
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Component B 
In our formula we must introduce the Geomagnetic field component in Moon’s direction. 


Referred to a system of axis North, East, Up: 

- The vector F is defined by Inclination and Declination. 

- Moon’s direction is defined by Azimuth and Elevation. 

For projecting Field F on the Moon's direction we must calculate the angle FM between these two vectors. 
The formula becomes: 


cosFM=cosl*cosD*cosEL*cosAz+cosl*sinD*cosEl*sinAz-sinl*sinEl 


from which B = F * cosFM 


TEC (Total Electron Content) 


TEC is a key parameter for describing Earth’s ionosphere. 
It is measured in TECU (TEC Units) = 10°° electrons/m? 
The number of TECUs represent the total number of electrons present in a cylinder of 1 m? of section, 


crossing the ionosphere in the wave’s direction. 


VTEC sources 
Scientific Stations post VTEC (Vertical TEC) diagrams for each day. Among these we found the Royal 
Observatory of Belgium (ROB), that publishes VTEC histograms with values every 15’, and keeps an archive 


for every day of the year. Measures are made by the observatory located in Dourbes (Belgium). 


Typical summer VTEC diagram: 


Vertical Total Electron Content in DRBS 


20 


VTEC (TECU) 
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Date: ] 2012-08-19 Time (UT) 
Station: @ppes 


Typical winter VTEC diagram: 


Vertical Total Etectron Content in DRBS 


VTEC (TECU) 
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Date: [2043-01-14 ‘| 2013—01—11 Time (UT) 
Station: © pres 


Image Source: ROB Royal Observatory Dourbes - Belgium 
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From Dourbes to other places 


Spatial variation of TEC: 
- Longitudinal variation: Global trend quite regular and correlated to the local solar time 


- Latitudinal variation: The TEC value, varies non-linearly from the poles to the equator (geomagnetic) 


2013-01-16 13:00:00 UT 
90° 180° 


Total Electron Content (TEC) 
-180" -90" 0" 


"180 oo 90° 180" 
lonospheric Range Error (L1) /m 


0.00 1.62 3.24 4.86 6.48 8.10 9.72 11.34 12,96 14.58 16.20 
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TEC / TECU 
This curve represents the latitudinal variation: 
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For our purposes we have simulated it with the algoritm: 
TECU variation = 0,02*LAT*2-2,5*LAT+95 

So for TEC value calculation: 

1 — Dourbes VTEC at same time 


2 — Magnetic latitude of station 
3 — With algoritm representing the curve we find the correction to be applied to Dourbes VTEC 


Slab thickness 

Data on the Dourbes site give VTEC and slab 
thickness. These represent the transformation of the 
real ionosphere, with changing electron content with 
height and total thickness over 1000 km, in an 
equivalent ionosphere constituted by a uniform density 
layer with known TECU and thickness. Going through 
this slab gives the same effects as going trough the real 
ionosphere. The advantage is that two numbers define 
this ionosphere, without the need for integration of a 
complex graph. 


4.0€+05 6.0E+0S B.OE*05 
N (om-3) 
Electron Density 


Fig. Schematic electron density height profile 


Oblique passage (Slant TEC Source: Journal of Atmospheric and Solar-Terrestrial Physics 


Ince the slab is crossed obliquely, the number of 


encountered electrons increases. 

So TEC in our formula becomes TEC = STEC = Ka*VTEC, with Ka given by: 

With Earth radius=6367 km, and lonosphere beginning at 100 km height, with h=Slab Thickness 
Ka =(SQR((64672+h)2-(6367*coSEI)2)-SQR(64672-(6367*cosEI)2))/h 


This is easily found with Pitagora’s theorem from 


Earth Path 
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Elev ation (deg) 


10&+08 


13 


1.2E+08 


Final formula 


We have now the data for the complete formula: 
@=k*(F*cosFM)*(VTEC*corr*Ka)/f? 
For f=144 MHz, k/f*=1,14 with F in Gauss. 


Wave plane rotation is controlled by these variables: 
— Angle between Geomagnetic. field and Moon direction (cosFM ranges from 1 to -1) 
— TEC (constant or changing slowly, 100% to 30%) 
— Moon elevation (oblique passage Ka from 1 to 3,7) 


Collecting on-the-air data 
We needed many sources of data, geographically spread. 


A big help came from René PE1L, who collected for us all the data published on line in LiveCQ in a file, 
accessible to us. 

We sorted them in an Excel sheet by date, spotter and spotted. 

Example: 18/08/2012 — DGOOPK — PE1L, data, pol and level eae 


[2 |date (Jute sig (pol call Elsgoner(= 
(1927 | 18082012) 802 -23) 130|/DGOOPK_ PEL 
[1923] 18082012, 804) -25 428|/DG00PK | PEIL 
[1937] 18082012, 808 | -23| 4148|DGOOPK |PE1L 
(2036) 18082012; 943) 21) 173/DG00PK —|PEIL 
(2039) 18082012) 945, -21| B\DGOOPK ss PEIL 
(2045) 18082012) 953) -23| 64|DGOOPK |PEIL 
(2050) 18082012) 955) -22 61/DGOORK |PE1L 
(2055) 18082012 957) -23) 174, DGDOPK |PE1L 
[2057] 18082012; 969, -21| SS BBIDGOOPK —PEIL_ 
(2060) 18082012/ 1001, -22| 58, DGOOPK PEL 
(2084] 18082012, 1016) “24, 82|DGOOPK |PE1L 
(2086) 18082012, 1018, 20) B3|DGOOPK PEL 
(2090) 18082012; 1020-21) BB|DGOORK —PEIL 
(2094) 18082012) 1022) -20) 85|DGOOPK |PE1L 
(2fon) 1a082012) 1024, -25|117|DG00PK —/PEIL_ 
[2105] 1g082012/ 1026-20, B7|DGOOPK —/PEIL 
(2435) 18082012; 1100-25] 134/DGOOPK —_|PEIL 
(2ta7| 18082012) 1102/22) 138,DG00PK | PEIL 
[2142] 18082012) 1104) -23) 46 DGOOPK |PE(L 
(21a5] 18082012) 1112/23) 153/DG00PK —_|PEIL 
(2148) 18082012/ 1114) -22/t61/DGOORK REIL 
(2251 18082012) 1330, -24) 66/DGOOPK |PE1L 
(2ane] 1g082012/ 1608) -23| 1 1|DGOOPK /PEIL_ 
(2412) 18082012 1510) “24, 17 | DGOOPK |PE1L. 
(2415) 18082012 1512) -25| 17/DGOOPK ——|PE1L” 
(2gat| 18082012/ 1514) -23| SS 15DGOOPK —_PEIL 
(2448) 18082012, 1520) -22| 8/DGOOPK |PEIL 
2456) 18082012 1622) -22| 4|DGOOPK |PEIL 


So we could examine a great variety of cases, and start to identify tendencies. 
Note: The polarization measured by MAP65 is the algebraic sum of spatial offset and of two ionospheric 


crossings: the Faraday rotation of the up going transmitted wave and that of the returning echo. 
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Amount of rotation in an interval 


As a first check on our formulas, we compiled an Excel sheet for calculating what happens in a time interval, 


generally not greater of an hour. An example: 


es a 0,070 0.48387 108,78 72.44 


L 3 0,067 0,483387 160,75 75,35) 
[71s a a ee 
| 12 | wKtuwe | sus3svr | bao | oS 
"bate [ure _[Tocattine [ Becimat_[AZQ)___[ETaeeeetol ¥TE ¥eEc-DRBS Tree con 
FAugust 3, 2012| 


(20 Faugust3,2012F 008 o:20 24 
CS) See ee a a 24.38 


We checked many cases with good congruence between computed and observed rotations. 


Global common moon trend 


Having now confidence in the basic correctness of formula and correction coefficients, we proceeded to build 


a new Excel sheet, covering the entire common-moon period. 
Our Excel sheet: 


12: 4 = 2 35 ry 08! f 2, $912 eZ 0 — 51 68, 2 68, 4 


ee ene 05543] S55 


13,95) 0564] 2311] 


ea) a 
a 


5 12.30 — ec 
46 


56 


307 16.17 0829 eT 0) p08 08 Zz! 4 [es : 


f -0) 2 79 4 4 

13.30 13.35) 163) 22,3 185 2.28 13, 68) i -0,6956 563.4 9,83 73,4 73,4 

=e =" a = 187 au . -0,7319 432.6) 8,60 83,7 83,7, 
+ 4 30 
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Partial checks were possible using the LiveCQ decoded periods. 

An example for SP4MPB spotted by PASFPQ on 16/12/2012. All data were computed for 30’ intervals. The 
choice of these two stations depended on having a LiveCQ decode running continuously for 44’. 

And the graphs we obtain from it (with superimposed the decoded graph): 


‘Vertical Total Election Content In ORBS 


SP4MPB was active from 13.58 to 14.42 utc. 


In this phase, TEC had a quick decrease, followed by a brief increase pre sunset, then decreasing from 
sunset to night. Calculated and real trend are coherent 


A second example is I2FAK calling CQ in ARRL EME contest, on 1/12/2012. 
Station is 828 SSE of spotter PASFPQ. 

Decoded pol showed 90° transitions when he changed tx form H to V. 
Graph corrected for this is the right one. 


Dec. Pol. Dec. corr. -90° 


22agetnh § 8 BB 
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Xpol antenna: Tx H and V 


Night conditions, 

with increasing 

Moon elevation 
VTEC 


18h 19h 20h 21h 22h 23h Oh 


Night TEC was practically constant, so pol change derived mainly from increasing Moon elevation. 
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Polarization trends 

When reception is difficult due to unfavorable Faraday rotation, one often wonders if this will change quickly. 
We decided to explore if tendencies could be found and used as a guide. 

So we created a set of 8 computations for a full moon pass of IK1UWL station with stations in 8 different 
directions, N, NE, E, SE, S,SW, W, NW, all for the same moonpass, on Dec 19" 2012, in which the Moon 
was visible from 11.00 to 23.00 UTC. In this diagram the graphs are disposed in the outer ring. 
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In order to separate the eventual tendencies from changes in the ionosphere (which occur over a full moon 


pass, but are negligible in a short period), we computed also the rotations with constant VTEC and slab 
thickness, i.e. for an invariable ionosphere. 


The relative graphs are disposed in the inner circle. 
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¢ Spotter IK1UWL, band 144 MHz, on Dec 19, 2012 — Moon 11.00 — 23.00 UTC 
« All graphs computed for stations in a rose of directions 
@=k*(F*cosFM)*(VTEC*corr*Ka)/f2 
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The first noticeable fact is the similarity at the beginning and at the end of the moon pass, whilst ionospheric 
changes give different evolutions in the central part of the moon pass. 

Evidently, at moonrise and moonset the rapid change of inclination (which influences both cosFM and Ka) 
causes quick variations of rotation. 

The overall shape of Eastward and Westward graphs is similar to a sinusoid, with slower variations in the 
central part. Specifically, at MR, cosFM variation of the station having MR, has a dominant effect, and is the 
cause, in the first hour approximately, of the pol rotation. The same effect can be seen at MS. 

For Northern stations the factor Ka (change of Slant TEC) is the dominat factor at MR and MS. 

For Southern stations the dominant effect, with IK1UWL spotting, is more tied to cosFM (angle between 
earth’s Geomagnetic field and Moon direction) due to the field’s condition above IK1UWL station, but also Ka 
show a similar effect. 

Generalizing, during hours in which TEC changes little (during the day and during the night), Faraday 
rotation is affected mainly by angle changes (Moon elevation and Moon direction respect Geomagnetic field). 


Changes tied to TEC occur mainly during sunrise and sunset. 
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Conclusions 


bd QSB of JT65 decodes: Is caused by focusing or defocusing of our beam going through the waves of 
the windy ionosphere. 

* Faraday rotation: There are three phases in a Moon pass: 
1 - In the first hours after Moon rise the rate of change of polarization is high. Causes: 
a) — change of angle FM between Moon direction and magnetic field b) — change in length of 
ionospheric crossing (slant coeff. Ka) 
2— Inthe central part of Moon pass changes in angle FM and coeff. Ka balance each other, so 
polarization changes depend mainly from ionospheric evolution (of Total Electron Content) 
3 — Inthe last hours before Moon set the rate of change of polarization is high for the same causes of 
phase 1 
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